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We consider ofl-diagonal correlation functions of impenetrable bosons on a lattice. By using the 
1— ( ' Jordan- Wigner transformation the one-body density matrix is represented as (Toeplitz) determinant 

, ^ , of a matrix of fermionic Green functions. Using the replica method we calculate exactly the full long- 

' ' range asymptotic behaviour of the one-body density matrix. We discuss how subleading oscillating 

' terms, originating from short-range correlations give rise to interesting features in the momentum 

, distribution. 
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I. INTRODUCTION 



Correlation properties of one-dimensional bosons with short-range interactions is a long standing problem of statis- 
tical physics of integrable systems Despite the exact integrability, the correlation functions are hard to calculate 
^ . due to a complicated form of the Bethe Ansatz wavefunctions and corresponding expressions for the matrix elements. 
Recently there has been a renewal of interest in these systems due to their relevance for current experiments with 
cold atoms. Experimental progress made it possible to create strongly correlated bosonic states in the gas phase 0,0 
^ and in an optical lattice j^i5|,j6|, and approach the Girardeau- Tonks regime of impenetrable bosons. The physics of 
\ bosons in this regime is similar in many respects to that of free fermions 0, since the strong repulsive interactions 
P3 ' effectively play a role of the Pauli principle. There is, however, an important difference due to bosonic quantum 
O " statistics of particles. The statistics reveals itself in off-diagonal correlations of bosons, such as the one-body density 
I ^ I . matrix or momentum distribution. The calculation of these observables in the Girardeau- Tonks limit from the first 

priciples represents an important question of statistical physics. 
CN| , Impenetrable bosons in one dimension is an exactly soluble model both in the continuum case and on a lattice. In 
^ ■ the latter case it is equivalent to the quantum XX spin chain: the two states of spin 1/2 on a given site correspond 
to the presence/absence of a boson. The correspondence with fermions was put forward in the early work by Lieb, 
Schultz and Mattis 8], where this system was mapped onto a system of free fermions on a lattice by the Jordan- 
Wigner transformation. Using this transformation it was possible to represent transverse spin correlation functions 
. (corresponding to the off-diagonal correlations of bosonic operators) in the form of determinant of a matrix with a 
\^ ' size given by the correlation distance. For a long time only an approximate behaviour of these objects was available. 
: Schultz has shown that at zero temperature correlation functions for large distances undergo a power law-decay. 

The explicit form of this power-law decay has been found by Efetov and Larkin (IQj and has been recovered later by 
\ Vaidya and Tracy using a method of mapping the quantum one-dimensional XX chain onto the two-dimensional 
^ ' classical Ising model. This technique has been succesfully employed in the works of McCoy, Barouch and Abraham 
' \ [l^ for calculating finite-temperature correlations. Recent numerical studies 13] of impenetrable bosons on a lattice 
^ in the presence of a confining potential have shown a high degree of universality of the power law decay found in the 
Q uniform case. 

^ In this paper we reconsider the problem of calculating off-diagonal correlation functions of impenetrable bosons on a 

^ lattice. We start with the Jordan- Wigner transformation and represent the one-body density matrix as a determinant 
• w^ , of a matrix given by fermionic Green functions. Next, we rewrite the determinant in terms of an average over the 
' ensemble of unitary random matrices and calculate its asymptotic behaviour by employing the replica method. Since 
5h ' its first use to derive the pair distribution function for random matrices |l4j , the replica method has been successfully 
applied to strongly correlated systems such as Calogero-Sutherland models [Tsl [l6j and impenetrable bosons in the 
continuu m llTl . The essence of the replica method in this context is exactly the same as in the theory of disordered 
systems jig. It consists of modifying the quantities being averaged so that they depend on a parameter n which is 
assumed to be integer. Once the average is performed, the result is recovered by a suitable analytical continuation to 
non-integer values of n. The main advantage of this method is that it leads to the asymptotic non-perturbative result 
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for the bosonic one-body density matrix including the leading power-law term plus oscillating contributions reflecting 
the short-range physics of the problem and lattice effects. This result allows one to determine the main features of 
the momentum distribution of impenetrable bosons on a lattice. 

The paper is organised as follows: in Section^]we use the Jordan- Wigner transformation to obtain the determinantal 
representation for the one-body density matrix of impenetrable bosons. We discuss how the results of Vaidya and 
Tracy are obtained by considering the leading asymptotic behaviour of the determinant. In Sections IIIII and IIVI we 
present the replica method and derive our main result, the full asymptotic expansion of the one-body density matrix. 
Section contains conclusions and prospects. The mathematical details are given in Appendices. 



II. ONE-BODY DENSITY MATRIX AS A TOEPLITZ DETERMINANT 

Let Urn, be creation and annihilation operators of bosons on an infinite lattice with sites labeled by the index 
771. In the limit of strong repulsive on-site interactions the Hilbert space is projected onto a subspace of states with 
at most one particle in each site. Using the Jordan- Wigner transformation 

0<l<m Q<l<m 

the bosonic operators are related to fermionic creation and annihilation operators, c]„, c,„ obeying standard anti- 
commutation relations 

{Cm, cj} = Sm.l] {Cm, q} = {cj„, c[} = 0. (2) 

The one-body density matrix is defined by the following ground-state expectation value: 



9i{R) = 7:{iam+B. + al^^ji){a„,+al^)) 



= i((c™+H + cJ„+^) n (-ir^' (c„-ct„)>, (3) 

m<l<m-\-R 

where we have used the relation 

(-ir?-'=(Q-ct)(Q+ct) (4) 

and the fermionic commutation relations 1^)1. 

For _R = the one-body density matrix is related to the mean number of particles per site (filling factor) v as 

51 (0) = (c^Cm) ^\=''^^' (5) 

while for R > 1 we use the relation Q for every factor in the product in Eq. ^ and represent the one-body density 
matrix as an expectation value of 2R operators: 

9l{R) — 2 (^™+-R-^m+-R-l^»"+-R-l ■ • ■ B,n+lAm+iBm), (6) 

where we have defined Ai = ci + cj, Bi = c; — c|. We now use the Wick theorem and take into account that free 
fermionic expectation values are given by 

{Am+lB,n} = Gl 

{Am+lAm) — {Brn+lBm) — 0, 

where the Green function of free fermions is 

Gl = ((c™+,+c|„+,)(c„-cL)) = y ^ (n, -i- n_, - 1) 
2 sm-Kvl 

= ; 5i,o. 7 
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Here Uq = n_q is the ground state momentum distribution of free fermions in a lattice with fiUing factor so that 



n„ + ri_ 



1, 

-1, 



\q\ < ni/ 
\q\ > my 



(8) 



Then the one-body density matrix © is represented as a determinant of R x R matrix 
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(9) 



The determinant ^ has identical elements along each of its diagonals and thus belongs to the class of Toeplitz 
determinants: 



=det[g. 



j-fcJi,/c=i,...,i?,; 



91 = ^ f g{q)e'''dq, 



(10) 



where g{q) is called generating function (symbol). It is easily shown from (0 that in our case the generating function 
is 



5((7) = e*«(ng + n_, -1) 



-iq 



1 



-iq 



1 giTTi^gjg 



(11) 



It has two jump discontinuities at Fermi points q — dnriy. The large size asymptotic behavior of Toeplitz determinants 
with such singular generating functions was a subject of extensive studies in mathematics (for overview and references 
see 01). It has been shown for a certain class of singular generating functions that the large size asymptotics of 
Toeplitz determinants is correctly described by the Fisher-Hartwig conjecture [l9|. In the case of generating function 
with two jumps as in Eq. (|ll|l the analysis of the asymptotic behaviour of the Toeplitz determinant © is subtle 
and the Fisher-Hartwig asymptotic formula remains a conjecture. Being applied, it yields the following long-distance 
behaviour of the one-body density matrix: 



/ o\ Poo 

9i(R) = — 

TT 



smTTi^ 



R 



R 



00, 



(12) 



where poo = 7re^/^2~^/^A~^, and A = exp(l/12-C'(-l)) ~ 1.2824271 is Glaisher's constant related to the Rieman zeta 
function. Numerical computation of the determinant (0 for R 100 and various values of the filling factor shows that 
Eq. 1)12(1 provides an excellent estimate for the dominant smooth behaviour of the one-body density matrix at large 
distances. Substracting this dominant contribution from the exact numerical expression reveals interesting oscillating 
corrections. These corrections are sensitive to the specific value of the filling factor and reflect an interplay between 
short-distance interparticle correlations and lattice effects. The Fisher-Hartwig conjecture is unable to capture this 
behaviour. In the next Section we present the calculations of the asymptotic long-distance behaviour of the one- 
body density matrix, based on the replica method which has been recently applied to study correlation properties of 
integrable models. We shall see that this method is able to provide us with the dominant behaviour H12|l as well as 
with the oscillating terms. 



III. REPLICA CALCULATIONS OF THE ONE-BODY DENSITY MATRIX 

Below we describe the replica method for finding the large R asymptotics of the determinant Q. By standard 
manipulations the determinant in (PJ is transformed to the following i?-dimcnsional integral: 

^^(^)=2i?!|,(2^l^«(^")"'n^fe)"2\n ' (13) 

where 

AflXe^'^) = A(e'^Se'«=,...,e^9«) = det[e^('=-i)'J']fe,i=i,...,fl,= [] {e^'^" - e'"^) (14) 

l<k<l<R 
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is the Vandermonde determinant. The notation of average in (|13|l is justified by the fact that the measure of integration 
coincides with the distribution of eigenvalues of random unitary matrices drawn from Dyson's Circular Unitary 
Ensemble [2(tI |. Denoting exp(iqi) = zi and exp(i7ri/) — v and using the representation we arrive at an equation 



9iiR) 



I R 



(r 




{r 


- vzi\ 








- VZl) 



y-r |1 ~ VZl\ |1 - VZl\ 
. . 1 — VZl 1 — VZl 



(15) 



The replica evaluation of the average begins with modifying (replicating) the function being averaged in H15|) 
by taking the 2n-th power of absolute values. Then the original expression is recovered by taking n ^ 1/2. It is 
however assumed throughout the calculations that n is an integer (and 2n is even) so that many simplifications occur. 
Eventually, the limit n ^ 1/2 will be taken with care and we shall describe it in detail. 

For an integer n, a straightforward algebra shows that we are dealing with an average of a rational function: 



1 /-i^ |1 - wi/P" |1 



vzi\ 



\i=i 



VZl 



1 - VZl 



,,2nR 



'yr {V - Zlf^ [V ~ Zlf^ \ 

l\ 1-vzi I -VZl I ^' 



(16) 



The crucial fact is the existance of a a duality transformation relating the R - dimensional integral (jl6|l to a 2n — 2 = m 
- dimensional integral: 



Zn{R)^^ / d"'xAlix)'[[x,{l-Xc)[il-x,)v + vx,Y 

^*Jm Jo 



where 



Sm. — 



c=l 



r2(i + c)r(i? + m + i-c) 

r{R + TO + 2 + c) 



(17) 



(18) 



The proof of the duality transformation is presented in Appendix 1X1 

Note that in the dual representation (|17|l large distance R appears only as a parameter, which makes the dual 
representation an excellent starting point for the asymptotic expansion of gi{R). The main contribution to the 
integral (|17|l comes from the end-points of integration. Let p be a number of variables in the vicinity of a;_ = 
and p' = m — p he a. number of variables close to a;.;- = 1. We shift the integration variables and approximate the 
integrated high-degree polynomial as 



R{i~v-^y 



+ R{1 - v^) ' 
The integration measure in H17|l factorizes as 



[(1 - Xc)v + VXc] ' 

[(1 - Xp+d)v + vxp+d] 



1,- 



(19) 



d^x A^„(x)n-T^c(l-:^c) 



c=l 



1 



i?(l - 
1 

R{l-v^) 



d(p+1) 



p'{p' + l) 



p' 



(20) 



d=l 



Then, summing over all possibilities to distribute m variables among saddle points and performing the remaining 
integration over fc, £,'d obtain 
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where the integrals are performed using the Selberg integration formula 



nOC P p 

•^0 c=l c=l 



(22) 



Combining them with the constant (|18|l we observe that each term in the sum (|21|l is proportional to the total 
combinatorial factor: 



r^(m + 3) 

r(p + i)r(p' + i) 



e^l n 



T{R + m + 2 + c) 

r(i? + 771 + 1 - C) ' 



where 



nLir(rf)nr^7r(b) 

nr=ir(c) 



(23) 



(24) 



is a familiar combination from previous studies of replica theories 0, 0| . 

Now we are at a position to take the limit n —^ 1/2. It is important to notice that for even m — 2n — 2 there 
is always a central term p ~ p' ^ m/2 in the sum (|21|l which provides a smooth non-oscillatory contribution. It is 
characterized by a maximal degree of replica symmetry breaking, since it originates from the saddle point for which 
the number of components Xd ~ x- is equal to the number of components Xb 2^+. This term will be retained in the 
limit 71 — !■ 1/2 (ill ^ — 1) and this provides an operational definition of the correct analytic continuation. 

Making a change of summation index from pto77i/2 + fc = 7i— 1 + fcwe factorize the combinatorial coefficient (|24|) 



as F^+'' ^ AnDi"\ where 



A. = n 



r(c) 



J-l r(2n+l-c)' 

c— 1 ^ ^ 



= n 



r(?i + c) 



r(7i 



^^r(n + i-c) r(7i + i-c) 

This leads to an expression 



fc > 0; 



(n) _ 



(25) 

(26) 
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(27) 



Here we have defined 



2t!-2 



C„(i?) = (-l)2"-ir2(2n+l) n 



T{R + 2n + c) 
V{R + 271 - 1 - c) 



(28) 



and the summation over k in Eq. H27I) was extended to all integers. 

Now we take the limit ri ^ 1/2 term by term in the series (|27|l . The only complication appears with the overall 
numerical factors C„ and A„. The analytical continuation of these factors is explained in Appendix^ where we show 
that A\ = V2/3oo and Ci = 1/R. We also have r(l/2 + fc)r(l/2 - fc) = 7r(-l)'=. We therefore find that 



gi(i?) = Z.(i?) = ^ 

2 TT 



sm TTV 
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The first several terms in this asymptotic expansion of the one-body density matrix read 

1 

1 cos27r7^_R 9 cosAnvR 
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32768 sin^ R^ 



(29) 



(30) 



where the leading term coincides with the result H12|) . As is seen from Eq. H3U|I . for small fc the coefficients in the 
series (|29ll rapidly decrease with increasing fc. However, they will diverge for large fc, which shows that we are dealing 
here with asymptotic, rather than convergent series. 
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IV. PERTURBATION THEORY 

In Eq. (|30|l . only the leading contribution of each saddle point was taken into account. To obtain the subleading 
terms one has to construct the perturbation theory around each saddle point. Using the shift of variables (|19|l the 
integrand in (|17|) is expanded in series in l/R. As a result, using the definition of the average 

(31) 



[f{i)9{a) , = (/(c)) (5(0) {f(o)^rf d^^fmu^U^ce' 

\ I p,p' \ / p\ I p' \ I P ip Jq 

each term in the series (|21ll gets multiplied by 
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c=l d=l 

Up to terms of the order of it is sufficient to put p = ni/2 = —1/2, since the leading contribution of other 

saddle points is of the order of 1/R^ or higher. The function i^„i/2,m/2(-R) is then evaluated using the following 
averages (see chapter 17 in 20j) : 



The result then is 
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(34) 



which leads to the following asymptotic expression for the one-body density matrix: 



/ Pec 

gi{R) = 

TT 



smTTi^ 



R 



coi -Kv 1 cos2-KvR , 

^"^2^"^^i^?^~R^ + ---' ^^^^ 



Our result reproduces correctly the continuous limit. In this case it is useful to introduce a length scale, the lattice 
constant a, so that x — Ra is the distance and nv/a — kp is the Fermi wavenumber proportional to the density. 
Going to the continuous limit corresponds to having the filling factor v tending to zero while the product kpx = irvR 
is kept fixed. Dividing gi by the lattice constant to have a correct normalization, Eq. (I30|l becomes 

g,{x) = {l/a)g,{x/a) = (l - — ^ - J^S^ + ' • ') (^6) 

\kFx\^ V 32(fcFa;)2 8 {kFxy ) 



in agreement with the Vaidya and Tracy result |21||. The sign of the oscillating cosine term in Eq. (|36|l has been 
corrected as discussed in |17| . 



V. CONCLUSIONS AND PROSPECTS 



The result (|35|) provides an exact analyical expression for the long distance behaviour of the one-body density matrix 
for impenetrable bosons on an infinite lattice. It is valid for all values of the filling factor z^, in particular it reproduces 
correctly the continuous limit v ^ The structure of Eq. H35|) is in accordance with the hydrodynamic expansion 
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conjectured by Haldane [23 • Equations H25|) . and (|26|) (at n = 1/2) provide the exact values for the non- universal 
coeficients of the leading oscillating terms. 

The oscillatory terms in H35|l reflect the physics on short distance scales, of the order of the mean interparticle 
separation. They are analogous to the Friedel oscillations in the physics of fermions and provide yet another mani- 
festation of fermionization. The interplay between these short-distance correlations and lattice effects is interesting, 
for instance, in the case of half filling {1/ — 1/2), where the particles tend to occupy every second site. This results in 
the oscillations of the one-body density matrix with the period of the lattice. 

The momentum distribution cannot be reconstructed from the asymptotic expression (|35|) for the one-body density 
matrix, as it requires the knowledge of the short-distance behaviour. However, we can predict some of the important 
features. For sufficiently small momenta their distribution diverges as an inverse square root of the momentum. The 
finite size would result in the finite peak value of the momentum distribution proportional to the square root of the 
total number of particles. For momenta of the order of 2kF = ^-nv j a the derivative of the momentum distribution 
has a jump (the distribution itself has a cusp) resulting from the underlying Fermi surface of the fermionized bosons. 
Weaker singularities exist for higher multiples of 2fc_F. It is again interesting to see that for the case of half filling 
the jump in the derivative occurs just at the border of the Brillouin zone and therefore cannot be observed: the 
momentum distribution has no singularities. 

These features can be probed in current experiments with cold atoms. The most promising method is expected 
to be a combination of Bragg spectroscopy with the time of flight technique. In this case a Bragg pulse is applied 
after switching off the trap, when the density of the cloud is already sufficiently small so that the number of scattered 
atoms (proportinal to the momentum disitribution) can be detected with precision |23j| . This technique can allow the 
observation of the singularities such as cusps in the momentum distribution originating from the oscillating subleading 
terms in the coordinate correlation functions. 

There are several directions in which the present work can be extended. An immediate question concerns the 
finite value of the interactions or, in other words the case of the soft-core bosons. This case is similar to the case 
of the Heisenberg spin chain, for which there are recent results based on the Bethe Ansatz solution [24]. It is not 
clear at the moment how to obtain the explicit behaviour of the correlation functions away from the free fermionic 
point (equivalent to the hard-core bosons) where the Tocplitz determinant expression exists for the one-body density 
matrix. In the case of impenetrable bosons, our work can be extended to the domain of time-dependent correlation 
functions which contains information about elementary excitations and are related to the responce of the system to 
external perturbations. 
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APPENDIX A: DUALITY TRANSFORMATION 

The average (|16ll is symmetric under the interchange of v and so that using the fact that z^"^ = z/ it is possible 
to bring Eq. (fT?)|l to the form 

2n{R) = Hl^^[]z,(l-«z-;)2"-l(z;-Z;)2«-l^ 

T-r2n I R 2n \ 

= ..^ooJi-„.„^„ ^-^^^iXiXliu^-mvc-z^)), (Al) 

Ui,...,U2n-l-*V, U2n^0 \i = lC=l / 

where we have "split" 2(2n — 1) points in the product using variables Uc for c = 1, . . . , 2n — 1 and Vc for c = 2, . . . , 2n. 
We have also defined two additional variables vi and U2n in order to represent the factor zi in the product on the 
same footing. 
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Recalling the definition H13|l for the average and using an obvious relation 

R 2n 



,V2n,Zl, 



.zr) 



1=1 c=l 



A2„(wi, . . . ,i;2„)Ai^(zi, ...,zji) 
we rewrite Eq. p6|) as a many-body correlation function 
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(A3) 



of fermionic creation and annihilation operators ip^{sc), ipitc) in the ground state of R + 2n fermions. Up to normaliza- 
tion, the ground state wavefunction is given by the Vandermonde determinants: {z\0) cx Aii^2n{z). The expectation 
value of the fermionic operators is given by the Wick theorem in the form of 2n x 2n determinant: 



(0|V^(mi) . . . V^(u2«)^^(f2n) . . . i^ivi)\0) 



with the following 2-fermionic Green function: 
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fi+2n-l 

G{uv) = (0|V^(M)7/;(f)|0) = 

p=0 



1 — uv 



R+2n 



(A5) 



We now take the limits U2n — > 0, wi — > cx) using G{u2nv) — G{Q) = 1 and G{uvi) ~ ^. The factor 

R+2n-i common to the elements in the first column cancels with the corresponding term in the prefactor of the 
second line of Eq. (|A1(I . Taking out a factor from each row we arrive at the following representation: 
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1 1 ... 1 



(A6) 



where G{u,v) ~ — — v) and we have used the relation Am(wi, . . . , u™) ~ 

- nr=i ""■^A™(ui, . . . , 

We take the singular limit Uc u and Uc ^ u by a standard procedure, described, for example in Expanding in 
non-zero elements of the last row and first column leads to expressing Eq. I|A6|I as a determinant of (2n — 2) x (2n — 2) 
matrix of partial derivatives: 



dudyG dudlG 



dudl^-^G 
dldl^-^G 



dl'^'^dyG dt^'-^dlG ... dl''~^dl''^^G 



(A7) 



Now we use the following integral representation in which we recognize the expression for the hypergeometric 
function: 



dudyG{u, v) 



dx x(l — x)[(l — x)u + vx 
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r(i? + 2n-2) 
r(i? + 2n+l)r2(2) . , 

r(i? + 2n-2) r(4) Zn + 6,Z,^,s), 



(A8) 



9 



where s = 1 — uv. Splitting again the variables s —> Sc = ^ — UcV, for c = 1, . . . , 2n — 2 the determinant (jA7|) is 
represented once again as a ratio 



ZniR) = 



2r2(2n-i)n':7'r(c) 

Fo{s2) Fi{s2) . 



r(i? + 2n+l)r2(2) 
r(i? + 2n-2) r(4) 

• ^2n-3(si) 
^2n-3(s2) 



2n-2 •p-r2n-2 ff. 



c=l 



A2„-2(s) 



-fo('S2n-2) ^l(s2n~2) •■• ^2ri-3(s2n-2) 



(A9) 



where Fk{s) = [{s — l)ds]^F{—R — 2n + 3, 2; 4; s). In this function k derivatives ds can be commuted through the 
factors s — 1 to the right, producing additional terms which has at most k — 1 derivatives. These terms do not 
contribute to the determinant, since they already appear in the k — 1-th column. Next, we use the known properties 
of hypergeometric functions to show that 



Fk{s) = (s-l)^9M"^-2n + 3,2,4,s) = 
r(i? + 2n-2) r(2 + /c) r(4) 



(-fl-2n + 3)fc(2), 



(s - lfF{-R - 2n + 3 + /c, 2 + /c; 4 + fc; s) 



T{R + 2n-2-k) r(2) r(4 + fc) 
Substituting this expression into ljA9|l and taking the limit s 



(1 - .s)'"+2"-iF(i? + 2n + 1, 2; 4 + fc; z) 
s = 1 — uv we get 



Zn{R) = 



(-1) 



2n-l 



nfa 



r(2)r(i? + 2n + l) 



2 ^J-^ r(l + c)r(3 + c)r(i? + 2n - l - c) 
F(0;s) F{l-s) ... F(2n-3;s) 



(1-.) 



(»i-l)(fl+4n-2) 



dsF{2n - 3; s) 



(AlO) 



where F{p; s) = F{R + 2n + 1, 2; 4 + p; s). For the derivatives of the function Fp{s) we have 
{N + 2n+l)k{2)k 



d^F{p;s) 



(4 + p)fe 



-F{N + 2n + I + k,2 + k;A + p + k; s) 



iJ^±^!^±3jiBl(^l _ ,yR-2n+i+p-kp^_j^ -2n + 3+p,2+p;A + p + k;s) 



r{R. + 2n + l + k) r(4 + p) 
r{R + 2n + l) r(2)r(2+j5)' 



(1 - / dxa;P+i(l-x)'^+i(l-sa;) 



fl:+2n-3-p 



In this expression we used the properties of hypergeometric functions |26| along with its integral representation. 
Substituting this result into IjAlOp yields 



ZniR) 



(-1) 



2n-l 



n r2(i 



r(i? + 2n + c) 



(1 + c)r(i? + 2n - 1 - c 



fl(n-l) 



»1 2n-2 
/ dxi... dx2n-2 (1 - Xif{l - X2f ... (1 - X2„^2)'""' TT a;c(l 



(1 - sx2)^+2"-3 2:2(1-312)^+^""'* 



xl'~\l-SX2)'' 



(1 - sa;2„-2)^+^"~^ 2;2n-2(l - sa;2„-2) 



i?,+2n~4 



^2:^:2(1 -*^2„-2)^ 



(All) 

-A2„-2(a;). 



The first factor in the integrand is anti-symmetrized to yield Vandermonde determinant A2„-2(l — x) 
The integration can be performed column- wise. This enables one to extract another factor — A2„-2(2;). Finally, 
setting m = 2n — 2 and \ — s = and interchanging v and v we obtain Eq. p7|l . 
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APPENDIX B: ANALYTICAL CONTINUATION OF PRODUCTS A„ AND C„(ii) 

The main idea behind evaluating the products Y[c=i fc such as (|25|l and (|28|l for fractional or negative values of 
the product limit n is to turn the products into exponentials of sums of logarithms of each factor and then use the 
integral representation for the logarithm of Euler's gamma function. It is possible to perform the summation 
under the integral and obtain an integral representation where n enters only as a parameter. The resulting integral 
representation can be analytically continued to a desired value of rt. In the case of A„ the procedure is described in 
detail in Appendix B of Ref. We cite here the final result: 

Poo = ^1/2/^2 = ixe^/^2-^I^A-^. (Bl) 

For the function Cn{R) defined by Eq. H28I) we have 



l'^ (_l)2nSff(2,+ i) = E(l-r(i? + 2n + c)-lnr(i? + 2n-l-c)) 



2ri-2 



c=l 







dt 

T 







dt 
T 



c=l c=l 
■ ^^^g-(2«+l)t _ p-(4n-l)t p-(2n-l)t „-t 



e"'^' 5 + 2n(2n - 2)e" 

(1-e-*)' 



(B2) 

For n = 1/2 the first term in square brackets is exactly equal to e~^', and we obtain the desired analytical continuation 



lnC.(i?)=/ ^(e-«'-e-*)=lnl (B3) 







Thus Ci{R) = 1/R. 
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